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SUMMARY 


The  main  achievements  between  June  30,  1981  and  June  30,  1982  are: 

1.  A  program  for  mathematical  solution  of  a  multi-ccrnpartmental  model  for 
simulation  of  uptake,  distribution,  and  elimination  of  vapors  having  a  capacity 
limited  elimination  pathway  was  prepared  for  the  Apple  II  Plus  computer  and 
tested  by  simulating  a  variety  of  trichloroethylene  exposures.  The  model  accom¬ 
modates  up  to  7  carpartments,  with  optional  linear  elimination  in  each  ccrtpart- 
ment  and  an  additional  non-linear  elimination  in  one  compartment. 

2.  Using  rabbit  tracheas  we  shewed  that  vapors  of  water  soluble  substances 
only  (not  gases  or  other  vapors)  are  retained  on  the  walls  of  respiratory  airways 
during  inspiration  and  are  desorped  during  expiration.  Because  of  this,  pulmonary 
uptake  measured  from  the  difference  in  concentrations  of  inhaled  and  exhaled  air 
does  not  represent  systemic  uptake,  and  concentrations  measured  in  end  exhaled 

air  do  not  represent  alveolar  concentrations ,  thus  explaining  vhy  experimental  data 
differ  from  data  predicted  for  water  soluble  vapors. 

3.  We  developed  a  method  for  determining  intrinsic  organ  clearance  in  small 
animals,  based  on  the  assumption  that  during  steady  state,  concentrations  in  non- 
excretory  organs  are  equilibrated  with  arterial  blood  and  the  concentration  in  re¬ 
gional  venous  blood  is  equilibrated  with  the  concentration  in  the  appropriate  ex¬ 
cretory  organ.  Intrinsic  clearance  is  then  calculated  from  the  mass  balance  across 
the  organ. 

4.  We  showed  in  rats  exposed  to  non-effective  concentrations  of  halothane  that 
halothane  metabolism  is  a  capacity  limited  process  and  that  anaerobic  metabolic  path¬ 
way  is  enhanced  when  oxidative  pathway  reaches  saturation.  Halothane  oxidation  is 
also  inhibited  in  simultaneous  exposure  to  isoflurane  or  methylene  chloride.  Nitrous 
oxide  has  no  effect  on  halothane  metabolism. 
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INTRODUCTION 


The  overall  objective  of  the  project  is  to  design  econo¬ 
mical  and  informative  testing  of  subacute  and  chronic  toxicity  of 
new  volatile  substances.  The  specific  objectives  are:  1)  to  pre¬ 
pare  a  mathematical  model  for  simulation  of  uptake,  distribution, 
and  elimination  of  vapors  having  capacity-limited  clearance  2)  to 
obtain  experimental  data  supporting  the  model.  3)  to  study  the 
factors  affecting  nonlinearity  of  clearance  (concentration  depen¬ 
dence,  interference  of  inhalation  of  other  vapors) .  The  project 
was  planned  for  two  years. 

This  report  describes  research  conducted  in  our  laborato¬ 
ry  during  the  first  year  of  the  project  (June  30,  1981  and  June  30, 
1982).  During  the  first  year  our  efforts  concentrated  on  two  items: 
I)  the  preparation  of  mathematical  solutions  of  a  multi-compart- 
mental  simulation  model  describing  uptake,  distribution  and  elimi¬ 
nation  of  inhaled  vapors  having  one  capacity-limited  elimination 
pathway.  2)  the  finishing  of  the  manuscript  of  the  monograph  en¬ 
titled  "Modeling  of  Uptake,  Distribution,  and  Elimination  of  Vapors" 
to  be  published  in  two  volumes  by  CRC  press  in  January  of  1983. 

The  principal  investigator  is  the  editor  and  major  contributor. 

One  copy  of  the  manuscript  of  each  chapter  authored  by  the  principal 
investigator  and  of  the  chapter  authored  by  Dr.  J.  Vlach  were  sub¬ 
mitted  to  AFOSR  in  January  1982,  when  funds  for  the  second  year 
were  requested.  These  chapters  include  data  generated  by  AFOSR 
supported  projects,  acknowledgement  of  which  is  made  in  the  preface 


-3- 


of  this  book 


This  interim  scientific  report  is  divided  into  two  parts: 
Computer  simulation, and  animal  studies. 

COMPUTER  SIMULATION: 


The  basic  program  for  Apple  II  computer  for  solving  the 
non-linear  multi-compartment  model  was  developed  in  collaboration 
with  Dr.  Vlach.  This  specific  model  accommodates  up  to  seven  com¬ 
partments  with  optional  linear  elimination  in  each  compartment  and 
an  additional  non-linear  elimination  in  one  compartment.  Non-linear 
elimination  is  considered  a  saturable  process  described  by  an  equation 
similar  to  the  Michaelis-Menten  equation  for  enzymatic  reactions 
(refered  to  as  capacity-limited  elimination) .  The  print-out  of  the 
computer  program  is  on  the  next  pages. 

USER'S  MANUAL 

DATA  STATEMENTS:  The  information  on  compartments  is  given  by  data 
statements  2-36. 

Statement  2:  Indicates  number  of  compartments. 

Statements  4-10:  are  reserved  for  the  description  of  the  compartments 
Each  statement  gives  information  on  one  compartment.  The 
order  of  compartments  is  arbitrary, with  two  exceptions! 

1)  the  first  statement  (statement  4)  must  relate  to  the 
compartment  with  capacity-limited  elimination.  2)  the 
last  statement  must  relate  to  lung  compartment. 

Three  items  of  information  are  given  in  each  statement. 
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Program  Non-Linear 


Program  can  have  non-linearity  defined  by  Michaelis-Menten  equation 
in  the  first  compartment. 


*  3PM  nmmpcp  nr  noMPOPTMCNTR 

•m - mXw*  ~ 

£  t  M  f 

3  pen  g,ox,c  ; 

4  DRTp  15.2,0,05 

5  OPTh  11.4,0.4.5 
0  CPTP  e.0,0,31.5 
7  OPTp  9.5,9,58 

9  nfiTP  24,0,853 

9  DP7P  4.2,0,0120 

10  DCTR  8,0,11.1 

15  REH  INITIAL  CONDITIONS, ( TIS. CONC. /PORT. COEFF .  '■ 

10  OPTP  0,0 

17  OPTP  3,0 

18  0PTP  8,8 
13  OPTP  0 

2^  REM  I.IMPX,KH 
20  OPTP  325,8 

P£M  RTPP  T  K1  MTWJJT^P 

775 — rjRTH — if — - - 

35  RFM  REQUIRED  PRECISION 
38  OPTP  0.01 

200  REH  1 TER RT IONS  TERHINP7E0 

205  DIM  IK 7  >,UN< 7 >,U0L0<  7 >,F<' 75,« 7,7 5 

210  OIH  G<7>,GX<7>,C<7>,DX<7> 

215  PIN  CU< 2 >,CO< 2 >,XK 3 > 

220  INPUT  "  WPNT  PRINTER  ?  >' V/N>" ;R* 

IF  R*  =  »v"  THEN  PRINT  "PR41" 

230  PRINT  General  information  bn  simulation  study 

i35  INPUT  "  OGTE:  7rTp$ 

243  INPUT  "OBJECT I UE:  "iB* 

245  INPUT  "  EXPOSED  SUBJECT:  "  iC$ 

250  INPUT  "rnMPOUND:  "iS* 

255  PRINT  :  PRINT 

'700  PEBT5  P — Rea'cTIn~and  print-dUt”  the  parameters 
285  PRINT  "NUMBER  OF  COMPPRTMENTS*" ,N 
270  prtut  "ft.'  t  0X(I''  iv  I  >" 

275  FOR  I  =  1  TO  N 
280  REPO  G<  I  >,GX<  I  >,C(  I  '• 

235  PR  I  NT  G<  I  >  ,  GX<  I  -.C'  I  > 

290  KIPVT  T 

295  REM  "INITIAL  CONDITIONS 

T0O  PRTNT  "TVJTTTfil  COW  I TT  HNS" 

305  FOR  I  =  1TO  N  STEP  2 

310  IF  I  =  m  then  REPO  "OLD1"  In  8QTQ  320 

315  REPO  MOLD*!  I  >,UOLDn  +  n 

320  IF  !  *  N  TWEN  PRINT  MOLD'  I  ■  :  G0T0  735 

325  05 1 NT  UOLDC I > .UOLEK  I  -  1 > 

T"?n  NPVT  ’ 

335  REPO  UMPX.XM 

340  PP T  NT  "  l.iMOXa " I.IMflV  ,  "  Wa "  VM 

345  REPO  N 

359  ps^n”  ''rrcs  3I2E  ’N  •  4IN.  " 

*rc;? 


NON-LINEAR 


SS5  PPM  INITIAL !3E  UOLTSGES  v' ?  S| 

379  FOR  I  =  1  TO  H 

37*;  «xi>  =  i;gld<i> 

380  HEAT  T 

385  REM  INITIALIZE  INTEGRALS 

Teg-  "rr;p-T . =  'T TEH? - 

395  All  I  >  =  0. 

400  NEXT  I 

405  REN  SET  JACOB I BN  WITHOUT  NONLINEAR I TV 

410  WMl's  H  -  1  ' 

415  FOR  I  *  1*T0  M 

420  FOR  J  a  1  TO  N 

425  A<I,J)  «  0 

438  NEXT  J 

435  NEXT  I 

440  pop  I  a  l  TO  NH1 

445  A<T,I>  *  SCI)  +  GXCI)  +  C<  I  >  /’  H 

450  p<"  I  ,  N  "-  a  -  0.--  I  ", 

455  B<N.I>  »  -  G<I) 

4C0  s  p<n,n'>  +  e<  i  > 

4G5  NEXT  I 

479  A<N*N>  =s  p<N,N>  +  GCN>  +  C<N>  •-  H 
475  PEN  SET  INITIPLTIHE  AND  STEP  COUNTER 

4587  *‘5: - 

405  1ST  a  8 
.490  OOLO  as  0. 

495  ww  a  H  /  2. 


509 

PEH  CO-OLD  CURRENTS  POP  INTEGRALS 

=:0=: 

510 

515 

C0<1>  =  0. 

COC  2  >  =  0. 

REM  TIME  STEPPING  STARTS  MERE 

520 

525 

RbfNf  "Pr^O11 

INPUT  "WANT  PRINTER  ?  <V/M> 

•1 

jR* 

w 

535 

.540 

TP  s  TyT?  Ti4P^  pPTMT  ''PR&I  11 

PRINT 

INPUT  "SIMULATION  DURATION  IN  HIN  = 

It  Tf 

i  !  S. 

eu>=; 

fuPiiT  "EXIT'S1  !PE  CCNTENTPOf inN  =  ” 

"it 

TFiPDT  "  pRinY  ftrtER  steR  “ 

iNr^_ 

55T 

T  s  T  4.  W 

580 

IS*r  s  1ST  +  1 

ppm  FRPnjpT 

UN1"  I 

570 

FOP  I  as  1  TQ 

N 

«r?«r 

I  lhY  T  a  2.  *  1  U'  T  '•  -  UOLO<  I  > 

589 

NEXT  r 

585 

ppm  newton 

ITERATION  STARTS  HERE 

FOR  I  =  1  tq  nmi 
ill  s  CUNCI)  -  UN< N > >  *  SCI) 

FC  I  >  =  -  Ui  -  SXC  I  )  *  !.!}■!(  I )  -  r<  I  )  *  ( !  N/  T  ••.  -  fl.'TV,  /  W 

cym\  =  cyn)  -i-  ::i 

i  next  t 

I  FCN)  a  FCN)  -  6<N>  *  CLINCH)  -  E)  -  CCN)  *  .'UN-'N'-  -  :j/N'»)  /  w 
REH  con  NQNLIN,  COXTR T rnT T nx 
:  Ul  =  XH  +  UC  1 )  “ 

CX  =  UHhX  *  i A'  1  )  /  Ul 
i  FC  1  )  =  F<  1  -  CX 

L_REH  MODIFY  THE  ENTRV  OF  THE  JPCCSIQN 

HI  s  i.J  1  A  III  - - - — 

1  DN  =  UMRX  *  XH  /  Ul 
PC  1 , 1  )  =  PC  1,1)  -  OOLD  DN 
i  OOLD  =  DN 

■  -FEH  SOLUS  THE  NEWTON  CQUQT70K1 
G0SUB94£ —  - 

"■.4kfl  uBTHin  new  uoltbges 

FOR  I  =  1  TO  N 
UNCI)  =  UNCI)  +  DXC  I  ) 

NEXT  I 

__REM  CHECK  CONUERGENCE 
rZR  a  0. 

FOR  I  a  i  TO  N 
R  =  PBS  COXCD) 

IF  R  >  ER  THEN  ER  a  S 
NEXT  I 

}J  *  ||  PRINT  ’’REDUCE  STEPS  I ZE  OR  TNCREBSE  PRECISION" 

IF  ER  <  EPS  THEN  ' GOTO  7S0 
GOTO  595 

PEM  ITERPTTONS  TERMINBTED 
USB  EgjggQI  RgSULTs"" 

_  REH  CQLOJLQTE  CURRENTS  QND  TNTEGRQl  S 

FOR  I  a  l  TO  N  - - 

UQLOC I  )  =  UC I > 

UCI)  a  UNCI) 

NEXT  I 

CUC  t  )  a  ([S’  —  { !([  N  )  )  *  3f  }  j  ' 
cue  2)  =  CX 

Ul  a  < COC 1 >  *  cue t >  5  *  HH 
JE  ill  >  0  THEN  Xlf'D  =  yjcj',  +  iii 
IF  Ul  <  0  THEN  X I C  2  )  =  XI C  2^  -  n’l 
*'  *  t  =  f  CO*'  2  '•  +  C!JC  2  )  )  *  ww 
XIC3)  a  XIC3)  +  Ul 
ccc  i  ’■  a  cue  i  ) 
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NON-LINEAR  (4) 


.  I 


REM  ITERATIONS  tcph I NATEOjTRRNSRER* RESULTS 
TJT  tct  <  KiPts  Turs  GOTO  ScQ 
355  PRINT 
86®  1ST  »  ® 

365  V/.  -  T  +  .1 
870  PRINT  "TIHE«"#TV. 

375  PRINT  " CONCENTRAT I ON  UALUES" 

380  FOR  I  *  1  TO  N  STEP  2 
885  IF  I  =  N  THEN  PRINT  'X  N  > 

090  IF  I  <  N  THEN  PRINT  I  >,W  I  +  15 
395  NEXT  I 

900  PRINT  "RATES":  PRINT  "PULH. UPTAKE  ELIMINATION" 

905  PRINT  CUC1>,CIX25 

910  PRINT  "TOTAL  AMOUNT":  PRINT  "PULH. UPTAKE  WASH-OUT" 

915  PRINT  XI<1>,XK2> 

920  PRINT  "TOTAL  ELIMINATED  < NOWLIN 5" 

925  PRINT  XK35 

330  IF  T  >  =  TE  THEN  SOTO  520 
935  REM  GO  BACK  FOR  A  NEW  STEP 
■940  SOTS'  555“ - 

945  REM  LU  DECOMP.  AND  SUBSTITUT.  Substitute  to  solve  Newton  equation 
W  FGPT'  =T  TOHI  (call  in  690) 

955  KP1  =  K  +  1 
960  FOR  J  =  KP1  TO  N 
365  A(K,J5  =  ACK-J5  '  ACK,K> 

97®  FOR  I  ■  KP1  TO  N 

975  A(  I*J5  =  A<I,J5  -  A<  I/rO  »  A<K,J5 

980  NEXT  I 

985.-  NEXT  J 

990  _  MPXT  K 

395  0X'<  1  )  =  F<1)  /  fie  1 .  i  •. 

1008  FOR  I  =  2  TO  N 

1005  IM1  a  I  -  1 
1010  SUM  a  0. 

1015  FOR  J  a  1  TO  IM1 

1020  SUM  *  SUM  +  A<  I  j  J  a  DX*'J"* 

1025  NEXT  J 

A<I,I> 


1030  OXCn  =  <F(  I )  -  SUM  5 
1835  NEXT  I 
1040  FOR  K  *  1  TO  MM 1 
1045  j  s  N  -  «' 

1950  IR1  a  I  +  i 

1 055  SUM  a  OX*'  I 

1060  FOR  J  =  IPt  TO  u 

1065  SUM  a  SUM  -  OX*'  J  >  *  a,- 1,  n 

1070  NEXT  J 

1075  oxen  =  sum 

1080  NEXT  K 
1085  RETURN 
1030  END 
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The  first  number,  G,  indicates  the  FLOW  of  the  vapor 
to  the  compartment.  The  value  of  G  equals  perfusion  rate, 
F,  of  the  compartment  multiplied  by  blood-air  partition 
coefficient,  Xbl/air.  The  only  exception  is  the  last  com- 
partment,  for  which  G  equals  alveolar  ventilation,  V  ^v. 

The  second  number,  Gx,  indicates  INTRINSIC  CLEARANCE 
of  the  compartment.  The  second  number  is  zero  in  the  first 
compartment  (statement  4)  if  the  clearance  is  capacity  li¬ 
mited.  The  second  number  also  equals  zero  in  any  state¬ 
ment  describing  a  compartment  with  no  elimination  pathway. 

The  third  number,  C,  indicates  the  CAPACITY  of  the 
compartment  to  retain  the  vapor.  The  value  equals  volume 
of  compartment  V,  multiplied  by  an  appropriate  tissue-air 
partition  coefficient.  The  exception  is  the  value  for  the 
last  compartment,  which  is  determined  by  equation  (1) : 

C.  =  FRC  +  1/3  V.  .  .  +  V.  X.  ,  .  +  QX.  .  ,  .  (1) 

lung  tid  lung  lung/air  bl/air 


Where  FRC  =  functional  residual  capacity,  V  ^  =  tidal  vo- 
lume,  viung  =  volume  of  lung  parenchymal  tissue,  Q  =  cardiac 
output  and  X's  are  appropriate  partition  coefficients  at  37C. 

Statements  16-19:  are  reserved  for  the  definition  of  starting  condi¬ 
tions  ; that  is ; for  concentration  values  in  each  compartment 
at  the  beginning  of  simulation.  They  must  be  listed  in  the 
same  order  as  compartments  are  listed  in  statements  4-10. 

Each  statement  accommodates  for  starting  conditions  of  two 
compartments.  To  start  the  first  exposure,  the  values  of 
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starting  conditions  are  equal  to  zero  or  to  the  normal  level 


of  vapor  in  the  compartments.  For  desaturation,  or  for 
an  additional  exposure,  the  values  of  starting  conditions 
equal  to  concentrations  of  the  vapors  in  the  appropriate 
compartments  divided  by  appropriate  tissue-air  partition 
coefficients  (refered  to  as  concentration  values) .  The 
concentrations,  which  represent  starting  points  of  the  si¬ 
mulation,  can  be  obtained  experimentally  (for  example, 
concentrations  measured  in  tissues  of  animals  sacrificed 
at  the  end  of  exposure  define  the  starting  conditions  for 
simulation  of  desaturation) ,  or  can  be  obtained  by  simu¬ 
lation  of  previous  exposures  as  concentration  values  ob¬ 
tained  at  the  end  of  the  immediately  preceding  exposure 
(see  output  data) . 

Statement  26:  describes  the  nonlinear  elimination  in  the  compartment 

described  in  statement  4.  The  first  number,  V  ,  indicates 

max 

the  maximum  possible  clearance  in  the  compartment  (assu¬ 
ming  no  flow  restriction) .  The  second  number  in  statement 
26,  Km  ,  indicates  the  vapor  concentration  in  regional 
(compartmental)  venous  blood,  at  which  clearance  equals  one 

half  of  V  . 

max 

Statements  29  and  36:  define  the  accuracy  of  the  calculation.  In 
general,  the  smaller  the  numbers,  the  greater  the  accu¬ 
racy  and  the  longer  the  calculation  time.  Inaccuracy  ap¬ 
parent  mainly  in  output  data  for  the  first  minutes  of 
simulation  is  reduced  with  time.  If  the  chosen  numbers 
in  statements  29  and  36  are  inadequate,  the  display  on  the 


screen  will  demand  that  the  value  for  the  STEP  (statement 
29)  is  reduced  or  PRECISION  (statement  36)  is  increased. 

It  is  advisable  to  start  wit  .  two  short  simulations:  in  the 
first  simulation,  the  STEP  =  0.5  (statement  29)  and  PRE¬ 
CISION  (statement  36)  hundred  times  smaller  than  the  si¬ 
mulated  exposure  concentration  (statement  545) .  In  the  se¬ 
cond  simulation,  STEP  is  reduced  by  half.  If  the  output 
data  in  both  simulations  is  the  same,  the  values  used  in  the 
first  simulation  were  correct.  If  there  is  a  large  diffe¬ 
rence  between  the  two  outputs,  then  STEP  must  be  reduced. 

IMPUT  STATEMENTS:  The  imput  statements  235-250  describe  the  experiments 
(Date,  objective,  exposed  subject,  compound) .  Statement 
540  gives  a  choice  of  exposure  duration,  statement  545  gives 
a  choice  of  exposure  concentrations,  and  statements  220,525, 
and  550  control  the  printer. 

Statements  220  and  525:  give  option  to  print  the  data.  If  the  first 
call  "WANT  PRINTER?"  (statement  220)  is  answered  "YES" 

(press  "Y"),  then  the  information  given  in  statements  1-36 
will  be  displayed  and  printed  (see  example  on  page  12)  .  If 
the  answer  is  "NO"  (press  "N"),  then  the  data  will  be  only 
displayed.  Statement  525  activates  the  printer  in  a  similar 
way  as  statement  220,  giving  option  to  print  the  output  data. 

If  the  answer  is  "Y",  then  all  data  which  appear  on  the 
screen  are  printed. 

Statement  550:  this  statement  controls  when  the  data  are  displayed  and 

printed.  The  calculated  data  *re  available  for  display  or  prin¬ 
ting  for  such  short  simulation  period  as  defined  by  "STEP" 

(Statement  29) .  However,  less  frequent  display  of  data  is 
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EXAMPLE  OF  PRINT-OUT 
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usually  prefered.  (STEP  is  usually  a  fraction  of  a  minute. . 

When  "PRINTER  AFTER  STEP"  is  displayed,  a  one  integer  number 
must  be  pressed.  This  number  determines,  how  many  steps  before 
the  data  will  be  displayed  or  printed.  For  example,  the 
display  of  calculated  data  is  requested  after  each  simulated 
minute.  If  the  STEP  (statement  29)  equals  0.5, the  answer  is 
2  (=  1/0.5).  If  the  data  print  is  requested  only  after  each 
simulated  hour,  then  the  answer  is  120  (60/0.5).  If  the 
STEP  (statement  29)  equals  0.2,  then  the  answers  are  5  and 
300  respectively. 

Statement  540:  displays  "SIMULATION  DURATION",  asking  for  information 
on  duration  of  exposure  or  desaturation  which  should  be  si¬ 
mulated  (in  minutes) .  If  repetitive  exposures  are  simulated, 
then  the  simulation  time  is  counted  from  the  beginning  of  the 
first  exposure. 

Statement  545:  displays  "EXPOSURE  CONCENTRATION",  asking  for  information 
on  exposure  concentration  which  must  be  expressed  in  weight/ 
volume  units.  If  exposure  was  interupted,  answer  "0"  to  ini¬ 
tiate  the  calculation  of  desaturation  data. 

OUTPUT  STATEMENTS:  (870-925) :  print  the  following  information  for  each 
simulation  period  selected  in  statement  550  (see  example  of 
print-out  on  pg.  12) : 

a.  first  the  actual  simulation  time  "TIME"  is  displayed. 
The  time  in  minutes  represents  the  time  interval  which 
collapsed  from  the  beginning  of  the  first  simulated 
exposure. 

b.  CONCENTRATION  VALUES  for  compartments  are  next.  The 

-13- 


W' 


vapor  concentrations  in  compartments  are  not  given 
directly.  To  obtain  the  vapor  concentrations  in 
tissues,  the  displayed  concentration  values  must  be 
multiplied  by  the  appropriate  tissue-air  partition 
coefficients.  To  obtain  vapor  concentrations  in 
regional  venous  blood,  the  concentration  values  must 
be  multiplied  by  blood-air  partition  coefficients. 

The  concentration  values  are  displayed  (printed)  in  the 
same  order  as  the  compartmental  parameters  were  given  in 
statements  4-10.  That  is,  the  concentration  value  for  com¬ 
partment  with  capacity  limited  clearance  is  displayed  first 
(first  line  to  the  3eft) .  Next  is  printed  the  concentration 
value  for  the  compartment  defined  by  statement  5  (first  line 
to  the  right) .  In  the  second  line  are  concentration  values 
for  the  compartments  defined  by  statements  6  and  7,  and  so 
on;  the  last  printed  concentration  value  equals  the  concen¬ 
tration  in  the  alveolar  air.  If  this  value  is  multiplied 
by  blood-air  partition  coefficient , then  the  concentration  in 
arterial  blood  is  obtained.  If  the  last  printed  concentration 
value  is  multiplied  by  lung-air  partition  coefficient,  then 
vapor  concentrations  in  lung  parenchymal  tissue  is  obtained. 

The  concentration  values  can  be  also  used  to  calculate 
rates  by  which  the  vapor  is  transfered  from  the  lung  to  the 
compartment  or  vice  versa.  To  calculate  this  transportation 
rate,  the  concentration  value  for  the  appropriate  compartment 
is  substracted  from  the  last  displayed  concentration  value 
and  the  difference  is  multiplied  by  appropriate  G-value  (given 
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in  the  corresponding  statement  4-10) .  If  the  value  is  po- 
sitive,  it  indicates  tissue  uptake  rate,  if  the  value  is  ne¬ 
gative,  it  indicates  tissue  wash-out  rate. 

c.  pulmonary  UPTAKE  RATE  (left)  and  ELIMINATION  RATE 
in  the  first  compartment  (right)  are  given  next.  Pulmonary 
uptake  rate  is  printed  because  of  its  importance.  Other¬ 
wise,  it  can  be  easily  calculated  by  subtracting  the  last 
concentration  value  from  exposure  concentration  (given  in 
statement  545;  notice  that  it  can  be  zero  for  desaturation) 
and  by  multiplying  the  difference  by  alveolar  ventilation. 

If  the  PULMONARY  UPTAKE  RATE  has  positive  signs,  it  represents 
pulmonary  uptake.  If  it  has  a  negative  sign,  it  represents 
EXHALATION  RATE  or  desaturation  rate. 

All  elimination  rates  can  be  easily  calculated  by  multi¬ 
plying  the  concentration  values  of  the  compartment  by  the  G 

•v 

characterizing  the  compartment.  The  elimination  rate  for 
capacity  limited  clearance  is  displayed. 

Finally, the  output  gives  numerical  integrals  of  pulmonary 
uptake  rates  and  of  the  elimination  rate  in  the  first  compart¬ 
ment  (capacity  limited  clearance) .  The  integrals  represent 
total  amounts  retained  and  eliminated.  The  integral  for  pul¬ 
monary  uptake  is  split  into  two  values:  1)  if  uptake  rates 
are  positive,  then  the  integral  represents  total  pulmonary 
uptake  (on  the  left) .  2)  if  the  uptake  rate  is  negative, 

then  the  integral  represents  the  pulmonary  wash-out  (on  the 
right) .  This  means  that  during  the  saturation  period,  the 
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integral  on  the  left  increases  with  time  and  the  integral  on 
the  right  remains  constant.  On  the  other  hand,  during  desa¬ 
turation  period,  the  integral  on  the  right  increases  and  the 
integral  on  the  left  remains  constant. 

The  last  value  is  the  integral  of  capacity- limited  elimi¬ 
nation.  This  integral,  which  represents  total  amount  eli¬ 
minated  } alway s  rises.  If  pulmonary  wash-out  integral  and 
elimination  integral  are  subtracted  from  the  pulmonary  up¬ 
take  integral,  the  body  burden  is  obtained. 

COMMENTS 

1*  The  consistency  of  units.  It  is  essential  that  all  para¬ 
meters  be  expressed  in  the  same  units.  The  concentrations  must  be  expres¬ 
sed  in  weight/volume  units.  For  example,  if  time,  volumes, and  flows  are 
given  in  minutes,  liters  and  liters  per  minute  respectively,  and  ex¬ 
posure  concentrations  are  given  in  mg/l^then  in  output  statements  con¬ 
centration  values  are  in  mg/1,  rates  are  in  mg/min.  and  total  amounts  are 
in  mg.  For  small  animals,  min.,  ml.,  ml/min.  and  ug/ml  can  be  prefered. 

For  small  exposure  concentrations  ,ug/l  or  ng/ml  can  be  prefered.  K  - 

•  m 

values  must  be  in  the  same  units  as  concentrations;  V  -  value  must 

max 

be  in  the  same  units  as  rates, and  intrinsic  clearances  must  be  in  the 
same  units  as  flows. 

2.  Computation  time  on  Apple  II  Plus  is  rather  long.  It 
takes  about  10  hours  of  computation  to  simulate  a  24  hour  exposure. 

3.  Modeling  of  repetitive  exposures  or  exposures  with  varying 
exposure  concentrations  can  be  done  with  the  program  if  the  exposure  is 


presented  as  a  number  of  consecutive  exposures  to  constant  exposure  con¬ 
centrations  (When  exposure  is  interrupted,  the  exposure  concentration 
equals  zero.) 

This  is  achieved  using  the  following  steps: 

a.  The  data  statements  are  prepared. 

b.  When  "SIMULATION  DURATION"  is  displayed,  the  duration  of 
the  first  exposure  period  is  pressed. 

c.  When  "EXPOSURE  CONCENTRATION"  is  displayed,  the  exposure 
concentration  for  the  first  exposure  period  is  pressed. 

d.  The  program  runs,  and  the  output  is  displayed  and/or 
printed  in  intervals  pre-scheduled  by  statement  STEP. 

e.  The  program  returns  to  statements  calling  for  informa¬ 
tion  on  the  second  exposure  period.  SIMULATION  DURA¬ 
TION  now  equals'  the  sum  of  durations  of  the  first  and 
second  exposure,  and  the  EXPOSURE  CONCENTRATION  refers 
to  the  second  exposure  period. 

f .  .  Repeat  steps  b-e  for  each  exposure  period,  making  simu¬ 

lation  duration  equal  to  the  sum  of  exposure  durations. 

In  the  computer  print-out,  the  time  periods  and  integrals 
are  added,  so  that  simulation  time  (the  first  line  in  each  output) 
gives  the  time  from  the  beginning  of  the  first  exposure  (  =  sum  of 
exposure  durations) }  and  the  integrals  represent  total  amounts  retained, 
exhaled  or  eliminated  during  the  whole  simulation  time,*  that  is,  from 
the  beginning  of  the  first  exposure  until  the  time  indicated  in  the 
first  line  of  each  print-out.  The  amounts  retained,  exhaled  or  eli¬ 
minated  during  each  exposure  period,  are  calculated  as  the  difference 
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between  the  corresponding  integrals  at  the  end  of  the  two  consecutive 
exposure  periods. 


4.  Scheduling  of  printer.  Rational  scheduling  of  printer 
makes  the  data  easier  to  survey  and  economizes  on  computer  time  and 
printer  paper.  The  output  data  is  needed  usually  more  frequently,  when 
the  values  are  rising  or  declining  rapidly,  while  less  frequent  output 
data  is  needed  if  values  change  slowly.  To  arrange  for  such  a  rational 
scheduling,  the  exposure  is  simulated  in  more  than  one  run.  After  each 
run,  the  calculation  is  returning  to  statement  540,  inquiring  about 
"SIMULATION  DURATION",  "EXPOSURE  CONCENTRATION"  and  "PRINT  AFTER  STEP”. 
For  example,  it  is  desirable  to  obtain  print-outs  after  each  minute 

during  the  first  ten  minutes  of  exposure;  after  each  ten  minutes* 
between  10  and  60  minutes  of  exposure^nd  after  each  30  minutes  between 
one  hour  to  eight  hours.  Then,  the  exposure  is  interrupted  and  print¬ 
outs  are  needed  each  minute  during  additional  five  minutes  of  desa- 
^uration.  Such  exposure  i^jsimulated  in  four  runs.  Assuming  that 
exposure  concentration  equals  7  and  step  equals  0.2,  the  imput  state¬ 
ments  for  each  run  are  as  follows: 

Run  Is  Exposure  cone.  7  Simulation  duration  10  Print  after  step  5 

Run  2 s  Exposure  cone.  7  Simulation  duration  60  Print  after  step  50 

Run  3:  Exposure  cone.  7  Simulation  duration  480  Print  after  step  300 

Run  4:  Exposure  cone.  0  Simulation  duration  485  Print  after  step  5 

5*  Simulation  of  changes  in  physiological  parameters.  If 
any  parameter  described  in  input  data  (statements  1-36) ,  changes  during 
exposure,  then  simulation  must  be  interrupted  and  the  program  returned 
to  the  beginning  (press  "RESET").  Changes  in  parameters  must  be  made 


in  appropriate  data  statements.  To  secure  continuity  of  simulation, 
the  final  concentration  values  in  first  simulation  must  be  used  as 
initial  conditions  for  the  second  simulation  (indicated  in  statements 
16-19) .  An  example  of  such  simulation  is  exposure  during  which  the 
pulmonary  ventilation  and  cardiac  output  rise  temporarily  because  of 
strain  excercise. 

6.  Preparation  of  data.  On  page  20  is  a  format  of  the  pro¬ 
tocol  we  are  using  for  preparation  of  data  for  simulation.  The  pro¬ 
tocol  sheet  is  self-explanatory.  For  demonstration,  a  protocol  sheet 
for  simulation  on  page  12  is  shown  on  page  21. 

We  have  by  now  available  a  series  of  simulations  of  dif¬ 
ferent  trichloroethylene  exposures.  These  simulations,  compared  with 
experimental  data,  confirm  the  validity  of  the  model. 
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ANIMAL  STUDIES 


1.  Retention  of  vapors  in  respiratory  airways.  The  simu¬ 
lation  models  developed  in  our  laboratory  reliably  predicted  the  bio¬ 
availability  of  vapors  and  gases  of  organic  solvents,  but  failed  to 
predict  the  bioavailability  of  vapors  and  gases  of  such  substances  as 
acetone  and  alcohols.  For  such  substances  as  acetones  and  alcohols, it 
was  apparent  from  comparison  of  predicted  and  experimental  data,  that 
transfer  of  vapor  from  ambiental  air  in  arterial  blood  was  less  than 
predicted.  We  suggest  that  discrepancy  between  predicted  and  experi¬ 
mental  data  is  caused  by  adsorption  and  desorption  of  the  substances 
in  the  walls  of  respiratory  airways.  This  study  was  undertaken  to  eva¬ 
luate  this  adsorption  and  desorption. 

Method : 

New  Zealand  rabbits  (weighing  about  2.5  kg)  were  used  for 
the  study.  The  rabbits  were  sacrified,  the  tracheas  were  immediately 
removed,  and  the  ends  of  the  tracheas  were  attached  to  swagelock 
T-pieces,  as  shown  in  the  scheme. 


2 

Tracheas  (usually  about  4cm.  long  with  inner  surface  area  6  cm  )  were 


wrapped  in  a  soft,  wet  cloth  and  kept  at  37C.  Sampling  syringes  were 
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50  yl  Hamilton  gas-tight  syringes.  Gas  and  vapor  concentrations  used 
in  the  studies  were  5  yM/1,  25  uM/1,  and  125  yM/1.  The  studied  gases 
and  vapors  are  listed  in  table  3.  From  the  side  of  the  preparation 
where  the  larynx  side  of  the  trachea  was  connected  with  the  T-piece  ^ 
20  ml.  of  gas  mixture  (equal  to  tidal  volume  of  a  rabbit)  were  pushed 
rapidly  through  the  trachea  from  a  glass  syringe.  Immediately  after, 
the  gas  samples  (25  yl)  were  collected  simultaneously  on  both  ends  of 
the  preparation  and  analyzed  by  gas  chromatography.  The  retention  in 
trachea  was  calculated  from  the  difference  of  concentrations  in  inflow 
and  outflow  ends.  The  retentions  in  table  3  are  expressed  as  a  frac¬ 
tion  of  the  injected  concentration. 


c  -  c 
in  out 


The  trachea  was  washed  with  five  strokes  of  clean  air  from  the  other 
side  of  the  preparation,  and  the  gas  mixture  was  injected  again.  The 
procedure  was  repeated  ten  times.  The  mean  value  was  used  for  evalua¬ 
tion  of  retention  in  each  trachea.  The  means  in  Table  3  were  calcula¬ 
ted  from  retentions  obtained  with  n-number  of  tracheas.  For  compounds 
which  had  large  retention  in  trachea,  desorption  was  observed.  Imme¬ 
diately  after  20  ml  of  gas  mixture  was  passed  through  the  trachea,  20 
ml  of  clean  air  was  passed  through  the  trachea  from  the  opposite  direc¬ 
tion,  and  gas  samples  were  collected  simultaneously  on  both  ends  of  the 
preparation.  Concentrations  of  the  studied  gas  were  measured  and  de¬ 
sorption  was  evaluated  from  the  difference  of  concentrations  across 
the  trachea. 
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C^ut~C|n  Where  and  C^ut  relate  to  inflow  and 

D  =  -  outflow  concentrations  after  injecting 

C.  clean  air. 

in 

For  substances  with  small  retention,  D  was  smaller  than  0.001.  How¬ 
ever,  a  significant  fraction  of  original  concentration  of  ethanol 
(D  =  0.2)  and  acetone  (D  =  0.1)  remained  on  trachea  walls  washed  with 
a  tidal  volume  of  clean  air.  Traces  of  ethanol  and  acetone  were  found 
in  the  fifth  wash  of  the  tracheas  with  a  tidal  volume  of  clean  air. 

This  study  is  finished. 

Conclusions 

The  retention  of  vapors  on  the  surface  of  the  rabbit  trachea 
is  shown  in  table  3.  The  absolute  value  of  the  measured  retention  is 
questionable,  since  only  retention  in  the  trachea  was  measured.  In 
vivo,  the  buccal  absorption  and  adsorption  in  other  parts  of  respira¬ 
tory  airways  would  increase  the  retention.  However,  since  retentions 
shown  in  Table  3  were  determined  under  the  same  experimental  conditions, 
the  values  indicate  the  relative  retention  among  the  substances. 

Our  model  is  based  on  the  assumption  that  no  vapor  or  gas  is 
retained  in  respiratory  airways.  Under  this  assumption,  the  concentra¬ 
tion  in  end-exhaled  air  equals  the  concentration  in  alveolar  air,  and 
pulmonary  uptake,  measured  as  the  difference  of  vapor  concentration  in 
inhaled  and  exhaled  air  multiplied  by  minute  ventilation,  represents 
systemic  uptake.  However,  if  adsorption  and  desorption  in  respiratory 
airways  take  place,  then  concentrations  in  end-exhaled  air  collected 
during  exposure  is  larger  than  alveolar  concentration,  and  pulmonary 
uptake  does  not  necessarily  represent  systemic  uptake.  These  studies 
indicate  that  our  existing  model  is  not  suitable  for  simulation  of 
uptake,  distribution,  and  elimination  of  water-soluble  compounds,  and 
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TABLE  3 


Retention  of  Vapors  in  Rabbit  Trachea  In  vitro. 


Percentage  of  Inflow  Concentration  Retained 

Substance 

n 

mean 

SD 

SE 

2,2, 2-trifluoro-l-chloroethane 

6 

-0.1 

3.3 

0.9 

1, l-difluoro-2-chloroethylene 

5 

0.0 

3.1 

1.0 

Methylene  Chloride 

5 

2.5 

9.5 

3.0 

Halothane  (CF3  -  CHClBr) 

5 

6.7 

1.5 

0.5 

Freon- 12  (dichlorodifluoromethane) 

5 

6.9 

!  1.0 

0.3 

Trichloroethylene 

5 

8.3 

5.0 

1.6 

Toluene 

10 

10.6 

5.5 

1.2 

Ethylene  Oxide 

4 

13.9 

4.1 

1.4 

Styrene 

6 

17.6 

6.0 

1.7 

Ethyl  Acetate 

9 

17.7 

9.9 

2.4 

Acetone 

9 

21.9 

7.5 

1.8 

n-Amvl  alcohol 

5 

41.1 

4.1 

1.3 

Acetyl  Acetone 

6 

45.8 

4.9 

1.4 

1-Butanol 

7 

48.2 

11.4 

3.0 

1-Propanol 

5 

54.0 

4.6 

1.5 

Ethanol 

11 

58.2 

11.4 

2.4 

1,4-Dioxane 

3 

59.5 

4.1 

1.7 

Methanol 

5 

68.6 

2.3 

0.7 

Retention  was  calculated  from  the  difference  between  inflow  and  outflow  concentra¬ 
tions: 


c.  -  c 
in  out 


R  =  100 


% 


in 


n  equals  the  mitiber  of  experiments  performed. 
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requires  further  modification. 


2.  A  new  method  for  determination  of  intrinsic  clearance 
in  vivo  was  developed.  Calculation  of  clearance  is  based  on  the 
concentration  difference  of  vapor  in  arterial  and  regional  venous  blood. 
The  concentrations  are  derived  from  vapor  concentrations  in  brain  (or 
muscle)  and  in  the  excretory  organ  at  steady  state  by  dividing  the 
tissue  concentrations  by  the  appropriate  tissue/blood  partition  coef¬ 
ficients.  This  method  is  simpler  than  the  previously  suggested  deter¬ 
mination  from  arterial  blood,  and  is  now  used  for  the  determination  of 
intrinsic  clearance  required  by  our  simulation  model. 

3.  The  thermodynamic  theory  for  estimation  of  partition 
coefficients  was  tested.  Stern  and  Shiah  (Molecular  Pharmacology,  19, 

56,  1981)  published  a  new  correlation  between  partition  coefficients 
and  critical  values  of  substances.  Since  determination  of  partition 
coefficients  is  cumbersome,  we  attempted  to  use  thermodynamic  theory 
for  prediction  of  blood/air  and  tissue/air  partition  coefficients.  The 
attempt  was  unsuccessful,  because  the  differences  between  predicted 
and  measured  values  were  too  large  (up  to  100%) . 

4 .  Dependence  of  halothane  metabolites  formation  in  rats 
on  exposure  concentration  and  on  presence  of  other  vapors. 

Procedure:  Sprague-Dawley  female  rats  (about  200  g)  were  exposed  for 

three  hours  to  halothane  or  to  a  mxiture  of  halothane  and  some  other 
vapor  (methylene  chloride,  isoflurane  or  nitrous  oxide) .  At  the  end 
of  the  exposure,  rats  were  sacrified  by  decapitation,  and  tissues  were 
excised  and  analyzed  for  halothane  and  its  metabolites  (trif luoroacetic 
acid,  1, l-difluoro-2-chloroethylene  and  1, 1, 1-trif luoro-2-chloroethane) 
and  for  the  presence  of  other  vapors  used  in  the  exposure.  The  exposure 


V 


procedure  and  analytical  method  were  the  same  as  in  our  previous 
project  (Biomathematical  Modeling  of  Chronic  Toxicity,  Progress 
Report  to  AFOSR  September,  1979) .  The  results  are  presented  in 
Figures  1-6.  Each  bar  represents  mean  and  standard  error  of  the 
mean  for  five  rats. 

a.  Dependence  oof  metabolites  formation  in  exposure 
concentration  is  shown  in  figure  1.  The  range  of  exposure  concen¬ 
tration  varied  between  70  p.p.m.  -  3,200  p.p.m.  (recommended  TLV 
is  50  p.p.m.,  anesthetic  alveolar  concentration  is  6,500  p.p.m.). 

In  the  graph,  concentrations  in  tissues  and  in  air  are  expressed 
in  yM/kg  or  yM/1  respectively.  These  units  permit  comparison  of 
the  formation  of  metabolites  on  the  equi-molar  basis.  Formation  of 
both  volatile  metabolites  (trif luoro-monochloroethane  and  difluoro- 
monochloroethylene  are  gases)  increases  sharply  with  increasing  ex-  ' 
posure  concentration.  Concentration  of  the  major  metabolite,  tri- 
fluoroacetic  acid,  remains  unchanged.  The  nonlinear  dependence  of 
metabolites  formation  on  exposure  concentration  is  more  apparent  if 
ratios  of  tissue  concentration  to  exposure  concentration  are  plotted 
against  exposure  concentration  (Figure  2) .  If  formation  of  metabo¬ 
lites  is  a  first  order  process,  then  the  tissue-exposure  concentration 
ratios  would  be  independent  of  exposure  concentration.  The  graphs 
show,  however,  that  concentrations  of  metabolites  in  tissues  increase 
less  than  exposure  concentration,  and  halothane  concentration  in 
tissues  increase  more  than  exposure  concentrations.  The  conclusion  is 
that  metabolite  formation  is  a  saturable,  flow-limited  process.  At 
exposure  concentration  25mg/l  (0.32%),  the  liver  exposure  air  concentra- 
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FIGURE  1 

CONCENTRATIONS  OF  NAIQTHANE  METABOLITES  IN  LIVER  OF  RATS 
EXPOSED  FOR  3  HOURS  TO  DIFFERENT  CONCENTRATIONS  OF  HALOTHANE 
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FIGURE  2 

DEPENDENCE  OF  HALOTHANE  METABOLITES  FORMATION 
IN  LIVER  ON  EXPOSURE  CONCENTRATION  OF  HALOTHANE 
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tion  ratio  approaches  the  liver-air  partition  coefficient  ^iiver/air® 
6.2).  This  indicates  that  at  this  concentration  metabolic  clearance 
can  be  considered  for  flow-limited.  It  is  suggested  that  the  enzyma¬ 
tic  system  is  saturated  if  exposure  concentrations  exceed  0.03%  (300 
p.p.m.),  since  the  concentrations  of  trifluoroacetic  acid  were  the 
same  in  tissues  of  rats  exposed  to  0.03%  of  halothane  as  in  rats  exposed 
to  0.32%  of  halothane. 

After  a  three-hour  exposure,  trifluoroacetic  acid  concentra¬ 
tion  in  liver  is  always  greater  than  halothane  concentration,  but  con¬ 
centrations  of  volatile  metabolites  are  always  smaller,  than  halothane 
concentration.  (Figure  1) .  Trifluoroacetic  acid  and  halothane  were 
also  found  in  other  tissues  in  a  significant  amount,  but  the  two 
volatile  metabolites  were  not  detectable  in  other  tissues  by  our  method 
(sensitivity  of  the  method  is  0.5  yM/kg) ,  The  explanation  is  that  the 
volatile  metabolites  formed  in  liver  are  removed  from  blood  circulation 
by  effective  pulmonary  clearance,  thus  preventing  distribution  of  the 
two  gaseous  metabolites  in  the  body.  On  the  other  hand,  the  trifluoro¬ 
acetic  acid  is  not  removed  in  lung  and  is  distributed  by  arterial  blood 
through  the  whole  body.  Since  renal  clearance  is  inefficient,  trifluo¬ 
roacetic  acid  accumulates  in  the  body,  and  after  a  short  exposure,  con¬ 
centrations  of  trifluoroacetic  acid  in  tissues  exceed  halothane  concen¬ 
trations. 

b.  Effect  of  simultaneous  exposure  to  another  vapor  on  for¬ 
mation  of  halothane  metabolites.  Rats  were  exposed  to  a  mixture  of  two 
vapors  in  which  halothane  was  always  present  at  a  concentration  5  mg/1 
(625  p.p.m.).  The  mixtures  of  halothane  with  three  vapors  were  stu¬ 
died  so  far.  Each  vapor  was  inhaled  in  four  sub-anesthetic  concentra¬ 
tions.  The  data  is  summarized  in  figures  3-6.  The  mean  tissue-concen- 
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tration  during  each  exposure  to  a  mixture  of  vapors,  and  standard 
errors  of  the  means, are  shown  by  bars.  The  mean  concentrations  in 
tissues  of  control  rats  exposed  only  to  halothane ,and  standard 
errors  of  the  means, are  shown  by  the  horizontal  lines. 

lb.)  Nitrous  oxide  was  present  in  the  mixtures  in  con¬ 
centrations,  3,10,30  and  50%  respectively.  Oxygen 
was  added  to  the  mixture  so  that  oxygen  concentra¬ 
tion  in  the  mixture  was  20%.  The  presence  of  ni¬ 
trous  oxide  had  no  effect  on  metabolites  formation 
(Figure  3) . 

2b.)  Isoflurane,  a  fluorinated  anesthetic  (CF3* CHC1*0,CHF2)/ 
is  metabolized  to  «  very  small  extent  (less  than  1%). 
The  molecular  ratios  of  halothane  to  isoflurane  in  the 
inhaled  mixtures  were:  1:0.'24;  1:1.1;  1:2.7;  and  1:5.2. 
In  all  exposures,  the  ratio  of  halothane  concentrations 
in  liver  and  inhaled  air  equals  four.  This  indicates 
that  systemic  clearance  of  halothane  is  flow  limited 
and  unaffected  by  isoflurane.  (Liver-gas  partition  co¬ 
efficient  equals  6) .  However,  isoflurane  affects  the 
individual  pathways  of  halothane.  Trifluoroacetic  acid 
formation  is  suppressed  in  the  presence  of  isoflurane. 
The  suppresion  is  concentration  dependent.  The  concen¬ 
trations  of  volatile  metabolites  are  slightly  increased 
in  liver  of  rats  exposed  to  those  mixtures  in  which 
isoflurane  concentrations  are  greater  than  halothane 
concentrations.  Even  if  the  concentrations  of  volatile 
metabolites  in  liver  indicate  only  slight  enhancement 


of  formation  of  volatile  metabolites  in  liver  of  rats 
exposed  simultaneously  to  isoflurane,  the  enhancement 
in  fact  can  be  large,  since  the  efficient  pulmonary 
clearance  prevents  cumulation  of  volatile  "etabolites 
in  tissues.  It  is  concluded  that,  quantitatively,  iso* 
flurane  does  not  affect  the  systemic  cleaiance  of 
halothane  to  any  significant  extent.  However,  isoflu¬ 
rane,  poorly  metabolized  by  itself,  inhibits  the  ae¬ 
robic  metabolic  pathway  of  halothane,  the  end  product 
of  which  is  trif luoroacetic  acid.  Isoflurane  enhances 
the  anaerobic  pathway  which  results  in  the  formatio'n  of 
volatile  metabolites  (CF3  -  CH2C1  and  CF2  =  CHC1) . 

Since  CF3  -  CH2Cl  is  hepatotoxic  (Brown  B.R.,  et  al. 

Env.  Health  Perspectives,  21,  185,  1977),  it  is  possible 
that  isoflurane  abet  hepatotoxicity  of  halothane. 

3b.)  Methylene  chloride  also  inhibits  metabolism  of  halothane 
(Figures  5  and  6) .  Results  of  the  experiment  with 
methylene  chloride  are  pictured  in  figure  5.  The  molar 
ratios  of  halothane  to  methylene  chloride  in  exposure 
mixtures  are:  1:1.1;  1:1.9;  1:4.5;  1:6.8  and  1:15.8. 
Halothane  concentrations  in  liver  show  that  systemic 
clearance  of  halothane  is  either  unaffected  or 
slightly  enhanced  by  the  presence  of  methylene  chloride. 
The  only  exception  is  exposure  to  a  mixture  of 
5  mg/1  of  halothane  and  30  mg/1  of  methylene  chloride 
(1:15.8),  during  which  systemic  clearance  seems  to  be 
reduced.  Concentrations  of  methylene  chloride  in  liver 
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correlate  with  exposure  concentration.  The  ratio  of 
methylene  chloride  in  liver  and  in  inhaled  air  oscil¬ 
lates  around  the  liver-air  partition  coefficient 
(Xiiver/air  =  suggesting  that  metabolic  clearance 

of  methylene  chloride  in  the  studied  concentration  range 
was  not  flow  limited.  The  formation  of  trifluoroacetic 
acid  is  largely  inhibited  in  the  presence  of  methylene 
chloride.  The  inhibition  is  concentration  dependent. 
Concentrations  of  volatile  metabolites  in  liver  of  rats 
exposed  to  mixtures  of  halo thane  and  methylene  chloride 
vary  widely  around  concentrations  measured  in  liver  of 
rats  exposed  only  to  halothane,  but  the  variation  is  not 
concentration  dependent. 

We  attempted  to  perform  similar  studies  using  ten 
times  smaller  concentrations  of  halothane  (Figure  6) . 
Metabolism  of  halothane  at  such  exposure  concentration 
(62  p.p.m. )  appears  to  be  a  flow — limited  process  of 
first  order.  Accordingly,  the  ratio  of  halothane  con¬ 
centration  in  liver  to  exposure  concentration  was  very 
small  (0.8  compared  to  liver-air  partition  coefficients 
equal  to  6) .  In  the  presence  of  methylene  chloride, 
this  ratio  increased  to  4,  indicating  reduction  of 
systemic  clearance  of  halothane.  At  the  same  time, 
formation  of  trifluoroacetic  acid  was  reduced,  so  that 
concentration  of  trifluoroacetic  acid  was  below  the 
sensitivity  limit  of  our  method.  Dif luoro-monochloro- 
ethylene,  which  was  present  in  liver  of  rats  exposed 
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only  to  62  p.p.m.  of  halothane  in  concentrations 
smaller  than  detectable  by  our  method,  was  measu¬ 
rable  in  liver  of  rats  exposed  to  a  halothane- 
methylene  chloride  mixture.  Concentrations  of 
trifluoro-monochloroethane  were  the  same,  whether 
halothane  alone  or  in  mixture,  was  used  for  ex¬ 
posure. 

We  conclude  that  methylene  chloride  supresses 
systemic  clearance  of  halothane,  profoundly  inhibi¬ 
ting  the  aerobic  metabolic  pathway  (formation  of 
trifluoroacetic  acid)  but  not  significantly  affecting 
the  anaerobic  pathway  (formation  of  volatile  metabo¬ 
lites)  .  We  speculate  that  the  inhibiting  effect  is 
caused  by  methylene  chloride  and  halothane  competing 
for  microsomal  mixed  function  oxidase,  or  by  desacti- 
vation  of  cytochrome  P-450  by  carbon  monoxide  released 
by  biodegradation  of  methylene  chloride. 

This  study  is  not  yet  completed.  We  are  initiating  similar 
experiments  with  mixtures  of  halothane  and  trichloroethylene. 
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